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We report on the application of a feedback-controlled electromigration technique for the forma- 
tion of nanometer-sized gaps in mesoscopic gold wires and rings. The effect of current density and 
temperature, linked via Joule heating, on the resulting gap size is investigated. Experimentally, a 
good thermal coupling to the substrate turned out to be crucial to reach electrode spacings below 
10 nm and to avoid overall melting of the nanostructures. This finding is supported by numerical 
calculations of the current-density and temperature proflles for structure layouts subjected to elec- 
tromigration. The numerical method can be used for optimizing the layout so as to predetermine 
the location where electromigation leads to the formation of a gap. 

PACS numbers: 66.30.Qa, 73.63.RT, 44.10.+i, 02.70.Bf 
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I. INTRODUCTION 

Electromigration leads to the formation of voids and 
hillocks in metallic conductors by thermally assisted ion 
diffusion under the influence of an applied electric field^. 
For more than a decade^, this process has been employed 
to deliberately form gaps of the order of a few nanome- 
ters only in microstructures. These gaps can be used 
to incorporate molecules in a source-drain channel. The 
planar geometry facilitates the integration of gating elec- 
trodes that are essential for a systematic characterization 
of electronic transport through molecules. 

Initially, electromigration was performed at cryo- 
genic temperatures and reported to yield irregularly 
shaped slits with a narrowest separation a. From low- 
temperature current- volt age measurements, a ~ 1 nm 
was estimated^'^. However, metallic clusters often 
remained in the generated gap^^. By working at 
room temperature, the cleanliness of the gaps could 
be improved, while the simultaneous implementation 
of a feedback- control still allowed the generation of 
nanometer-spaced electrodes^^. We used a realization 
of this latter technique aimed at the control of the power 
dissipated while a gap is being formed. The procedure 
is similar to the one described in Ref. 10. Our aim is to 
extend the technique used so far for gaps in single wires 
to more complicated structures, e. g., to form gaps in the 
two arms of a ring connected by leads at opposite sides. 

To obtain a high yield of nanogaps, the experimental 
parameters have to be optimized with respect to two com- 
peting objectives. On the one hand, the temperature T 
and the current density j have to be raised sufficiently for 
material transport to occur on a reasonable time scale. 
On the other hand. Joule heating has to be limited to 
avoid melting of the metal structure, which would lead 
to droplet formation due to surface tension and to an 
unwanted increase in electrode spacing. 



In our experiments, we find a clear dependence of the 
resulting gap size on the thermal coupling of the electro- 
migrated structure to the cold substrate. We present a 
numerical model which leads to a plausible explanation 
of this finding. It maps the current density and the ex- 
cess temperature in the voltage-biased structure taking 
into account the heat flow to the substrate. The cal- 
culations can be used to optimize the sample geometry 
so as to predetermine the location where electromigra- 
tion leads to the formation of voids and the structure is 
finally going to break. 



II. EXPERIMENTAL OBSERVATIONS 

All samples used for this work were prepared by 
standard e-beam lithography on commercially available, 
slightly p-doped (boron, resistivity: 1 • • • 20 ^^cm) silicon 
substrates. Their surface was either covered by 400 nm 
thermally grown Si02 or merely by native silicon oxide. 
Here, we focus on samples made from 30 nm thick gold 
deposited by electron-beam evaporation without any ad- 
ditional adhesion layer. 

Electromigration was performed at room temperature 
and mostly at ambient pressure. Some experiments were 
carried out in the high vacuum (10~^Pa) of a scanning 
electron microscope (SEM) for in-situ monitoring of the 
electromigration process. The resulting gap size, how- 
ever, did not depend on the background pressure. The 
active feedback- control loop was realized by ramping a 
voltage across the sample and measuring its resistance 
in a four-point configuration. Upon detection of a cer- 
tain change of AR/R ~ 5 ... 10%, the voltage was set 
back to a safe value immediately by the control software 
and the voltage ramp was reinitialized. At an advanced 
stage of the process a sudden change of resistance by sev- 
eral orders of magnitude indicates rupture of the metallic 




Figure 1 . SEM images of gold wires after electromigration in 
high vacuum (10~^ mbar). The substrate consists of p-doped 
sihcon covered by (a) native SiO^ and (b) 400 nm thermaUy 
grown SiOcc, respectively. 



current path and the procedure is stopped. 

In this paper we want to address the issue of Joule 
heating during electromigration. Its effect manifested it- 
self mainly in two aspects of our experiment. On the one 
hand, SEM imaging after gap formation revealed that 
metallic structures on silicon wafers covered by a 400 nm 
thick thermally grown silicon oxide layer tend to show 
signs of melting. Consequently, the resulting electrode 
separation is considerably above 10 nm. A typical ex- 
ample is shown in Fig. 1 (b). Samples electromigrated 
on silicon covered by native oxide, however, display the 
intended formation of voids and hillocks leading to the 
opening of an irregular gap, whose expansion lies below 
the resolution of the SEM, i.e., a < 10 nm. An exam- 
ple is shown in Fig. 1 (a). On the other hand, as can 
be seen in the same figure, the samples were usually dis- 
rupted near their cathode end and not at the constriction 
prestructured for this purpose at the center of the wire. 

Breaking at the cathode side is a common observation 
in straight wires^'^^^^. It arises because electromigration 
is due to two mechanisms leading to atom transport in 
opposite directions. First, atom transport it is directly 
driven by the electrostatic force acting on the positively 
charged atom cores in the applied electric field. This 
mass transport is directed towards the cathode. Second, 
the cores are driven by the so called ''''electron wind^\ i. e., 
the momentum transfer from the conduction electrons to 
defects in inelastic scattering events. The net contribu- 
tion of these impacts is parallel to the electron current. 
In gold the motion due to the electron wind predomi- 
nates and thus material is effectively transported towards 
the anode. Since the amount of material flow depends 
strongly on temperature, electromigration is enhanced 
by Joule heating in those parts of the sample where the 
current density is highest and leads to depletion of ma- 
terial at the cathode side of narrow constrictions in the 
current path. 

The purpose of the constriction shown in the center 
of Fig. 1 (b) is to enhance the current density and cre- 
ate a hot spot at a predefined position. Electromigation 
was expected to lead to a break junction at this very lo- 
cation which, however, was not observed. Apparently, 
the conditions at the lead-nanowire-j unction have to be 
taken into account more carefully. In the next section 
we present a model to calculate the temperature profile 



for a given sample layout which takes into account the 
heat flow from the metal film to the leads as well as to 
the substrate. The latter part is of utmost importance 
as indicated by our experimental observations presented 
in Fig. 1. 



III. CALCULATION OF THE HEAT 
DISTRIBUTION 

The temperature in a voltage-biased sample is en- 
hanced due to Joule heating. We present a model to 
quantify this effect. It is applicable to metallic structures 
of constant thickness d, sheet resistance R[j , and thermal 
conductivity /c^. Heat is carried along the metal film 
to the leads and through an insulating barrier of thick- 
ness b and thermal conductivity /C5, where the condition 
kb <C km naturally holds. The structure is supported be- 
low the barrier by a thick substrate with thermal conduc- 
tivity kg. The substrate is assumed to be thermally well 
coupled at the back side to a heat reservoir at ambient 
temperature. It is heated from the top by the power flow- 
ing through the barrier which results from the dissipation 
in the metallic structure. The situation in the substrate 
is therefore intrinsically three dimensional. However, we 
can safely assume that the temperature in the metallic 
structure itself does not vary significantly perpendicular 
to the film plain (owing to k^ <C km)- The tempera- 
ture rise in the substrate depends strongly on the bar- 
rier thickness b and is significant for small b only. We 
postpone the discussion of thicker barriers to the end of 
this section and assume first that the temperature varies 
linearly across the insulating barrier, i.e., we consider 
a laminar heat flow across the insulator perpendicular 
to its interfaces with the metal film and the substrate. 
The numerics for the metallic film is reduced to a two- 
dimensional problem obeying the modified Poisson equa- 
tion: 



k^dV'T + q=^{T-T,). 



(1) 



Here T(r) is the temperature elevation relative to the 
base temperature of the unbiased sample, q = j^Ro the 
power density, j(r) the two-dimensional current density, 
and Ts (r) the temperature at the interface between of the 
substrate and the barrier. The left-hand side of Eq. (1) 
describes the heat balance in the metallic film, while the 
right-hand side is equal to the loss due to heat flow to the 
substrate. The latter depends implicitly on Tg(r) which 
itself obeys the three-dimensional Laplace equation (de- 
scribing the source-free heat flow in the thick substrate) 



V^T, = 



(2) 



subjected to mixed boundary conditions. At the interface 
with the insulating barrier we have 



bk. 



n 



(3) 



where Vn is the gradient perpendicular to the interface. 
At ah other boundaries we approximate Ts(r) by zero. 
Equations (1) and (2) form a set of mutuahy dependent 
relations. 

We find it convenient to solve this set of equations by 
a finite- difference method on a square lattice and dis- 
cretize the relevant parts of the metallic structure with 
a predefined cell size sufficiently small to capture all fine 
details. For this purpose we use an appropriate rectan- 
gular cutout of the e-beam lithography design file, which 
covers the area of highest current densities. For the ex- 
amples presented below, a cell size between 1 nm and 
5nm turned out to be reasonable, with several million 
cells within the cutout. The current and voltage leads (cf. 
horizontal and vertical structures, respectively, emanat- 
ing from the center in Fig. 2) have constant cross section 
and extend beyond the cutout. In addition the length 
of the leads included in the explicit calculation exceeds 
their width. These measures guarentee that the equipo- 
tential lines beyond the cutout generated by an electro- 
motive force between source and drain current leads run 
perpendicular to the latter. 

Before solving Eqs. (1) and (2) we have to calculate 
the current density j which enters the second term on 
the left-hand side of Eq. (1). j is directly related by 
j = W/Ru to the potential V ^ which in turn obeys the 
Laplace equation V^^ = with mixed boundary condi- 
tions. We solve the latter equation by successive overre- 
laxation (SOR)^^. For edges of the metallic structure in- 
ternal to the explicitly treated area as well as for sections 
of voltage probes with its border we have von-Neumann 
conditions, VnV = 0. At the section of the source and 
drain leads with the border we apply the actual bias Vmax 
and ymin, respectively. After solving the Laplace equa- 
tion subjected to these boundary conditions, the total 
current / = j ds\VnV\/Ru can be calculated, where the 
integral runs along an arbitrary line cutting the source- 
drain path. It is convenient to normalize V by IRu^ 
since V/IRu directly reflects the '"''number of squares''' 
between different points within the metallic structure, 
i.e., the length-to-width ratio of a strip which would ac- 
quire an equivalent voltage drop if biased by the same 
current. Figure 2 shows the result of the calculation for 
a layout typically used in our experiments. The inset of 
this figure displays the power density in the center of the 
structure. The notch in the middle does indeed lead to 
the intended peak of the power density at this position 
(and of the current density that drives electromigration) . 
However, a similar peak is found at the crossover from 
the funnel-shaped leads to the narrow wire as indicated 
by the four lightly colored spots at the rectangular cor- 
ners in the inset of Fig. 2. This is due to the well-known 
divergence of the electric field strength at corners^^. In 
our numerics, this divergence is limited only by the finite 
cell size used in the finite-difference method as corners 
are infinitely sharp in the layout. In practice, e-beam 
lithography leads to a minimal radius at corners set by 
the finite resolution of the process and the granularity of 
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Figure 2. Equipotentential lines for a test layout biased in 
horizontal direction. The taps pointing up- and downwards 
represent voltage probes. The voltage drop between consecu- 
tive lines is Al/ = O.lRu • I and the total drop between the 
outermost lines V — 12 ARu ■ I. The inner region where the 
potential lines merge to a continuous color scale corresponds 
to /SV = 5i?n • I- The inset displays the resulting power den- 
sity q in units of I^Ru/w^ fo^ the center of the main figure. 
Notice the peaking of the power density at concave corners. 



the crystalline metallic layers. This radius sets the scale 
of the actual current-density enhancements. It can thus 
be expected that the current density reaches the same 
order of magnitude at each concave corner. The precise 
value depends, however, on the opening angles and can 
be tuned by proper design. Consequently, rectangular 
corners as the ones visible in Fig. 2 should be avoided in 
the strongly heated parts of the structure, unless a start- 
ing point for electromigration is intended. We are going 
to comment further on this issue in the next section. 

Having q = (VF)^/i?n ^it hand, we solve Eqs. (1) and 
(2) simultaneously by an iterative procedure. Starting 
with T(r) = Ts(y) = we first relax Eq. (1) at constant 
Ts to get an updated estimate for T(r). The bound- 
ary conditions for the structure edges are again of von- 
Neumann type, Vn^ = 0. Special care has to be taken 
to treat the heat flow at the border of the cutout. Owing 
to the aforementioned homogeneous nature of the leads 
extending beyond this border, the heat-flow problem is 
essentially one dimensional here^^. In long leads the tem- 
perature along the wire direction approaches exponen- 
tially a steady-state value of 

Ta = Ts^IJRuh/{w''h), 

where Ii is the current in the specific lead of width w. 
In general, the substrate temperature Tg in this equation 
has to be determined self-consistently as it results from 
the heat flow caused by the last term. However, T^ is 
naturally small compared to more strongly heated zones 
below spots of high current density. We ignore small 
variations of Ts at the border of the cutout and take in 
the above formula for Tg a self-consistently calculated 
value. The functional dependence of the temperature 



along the leads is then given by 

T{x) = T, + (To - Ta) exp(-(x - x^)/\). (4) 

Here x is a spatial coordinate along the current direction 
entering the cutout at xq, Tq is the te mperature of the 
cross section at xq, and A = ^/btkm/kb is the length scale 
of relaxation towards the asymptotic T^. From Eq. (4) 
the boundary condition of Eq. (1) at the cutout frame 
can be deduced, 



V„T = {Ta - To)/\. 



(5) 



Eq. (5) completes the boundary- value problem defined 
in Eq. (1) and is used in a successive overrelaxation loop 
to update T(r). In the next step of the iterative pro- 
cedure we fix T(r) in Eq. (3) and update Ts{r) via Eq. 
(2). Here, successive overrelaxtion would fail since the 3d 
nature of heat relaxation in the substrate leads to an in- 
tractably large number of cells in the discretization step. 
The final goal of our calculation, however, is the deter- 
mination of T(r) defined in the 2d plane of the metallic 
structure. To this end we need to know Ts{r) at the in- 
terface between the substrate and the barrier only. Let 
G(r) be the solution of V^G = in the 3d half- space of 
z < obeying Dirichlet conditions G{r -^ oo) = and 
von-Neumann conditions at z = 0, VnG = S{r). Due to 
the linearity of the Laplace equation we can write down 
immediately the solution of Eq. (2) obeying the boundary 
condition of Eq. (3), 



Tsir) = I; /dr'G(r - r') (T(r') - r,(r')) . 



(6) 



The integral is restricted to the plane z = 0. Thus, the 
calculation of Tg(r) at the interface z = requires the 
knowledge of G{r) at the interface only as well. In this 
way Eq. (6) is reduced to a two-dimensional convolu- 
tion which can be solved numerically by discrete Fourier 
transforms (DFT) for which fast algorithms exist. Our 
problem is thus reduced to the determination of G(r). 
Details of this calculation are given in the appendix. 

Eq. (6) can now be used to update Ts{r). Note, how- 
ever, that Eq. (6) does define Tg only implicitly. There- 
fore a further relaxation loop is required. In each itera- 
tion of this loop the convolution on the right-hand side 
is calculated employing the convolution theorem by two 
consecutive DFT steps. The result is subsequently used 
to update Tg. In the course of the iteration, the changes 
in Tg get smaller and the loop is terminated after suf- 
ficient accuracy has been obtained. This ends the first 
round of our overall iterative procedure. In the follow- 
ing iterations we proceed by using Eq. (1) at constant 
Tg to update T and Eq. (2) at constant T to update Tg. 
The procedure converges to a self-consistent solution af- 
ter several iterations. 

So far we have developed a model applicable for thin 
barrier thickness b. From now on, we refer to this partic- 
ular approximation as Approx. I. In the next section we 



present results showing that indeed substrate tempera- 
ture Ts and T follow each other closely for small b and 
the a-priori approximation of a laminar heat flow per- 
pendicular to the interface between substrate and barrier 
is justified in this case. However, Tg is a monotonically 
decreasing function of b while at the same time T rises 
monotonically Moreover, at small b it is possible to heat 
the substrate rather locally by the large power density 
resulting from a notch as, e.g., the one in the center of 
the test structure displayed in Fig. 2. If 6 is increased 
the power is delivered to the substrate more evenly and 
the temperature distribution tends to smooth out (with 
a small enhancement only even at the center). The tem- 
perature difference T — Ts soon becomes so large that Tg 
on the right-hand side of Eq. (1) can safely be neglected. 
In this case, the vertical nature of the heat flow is lost 
and the problem of the heat distribution in the isolating 
barrier has to be treated more precisely. Fortunately, this 
can be done by the same method introduced so far. For 
thicker barriers we take the temperature at the interface 
of the substrate with the barrier as unelevated, the bar- 
rier now stays cold at its back side and takes thus the 
role formerly played by the substrate. The set of Eqs. 
(1), (2), and (3) now reads 



kmdW^T ^ q = h- 



AT 



V^T, = 



-f, .) 



where AT is the temperature drop across a slice of the 
barrier of infinitesimal thickness e. In the finite- difference 
treatment of the problem we set e equal to the discretiza- 
tion cell size. In Eq. (7) Tg = Tg(r) represents the 3d 
temperature profile in the barrier. As such, it is defined 
at the interface of the barrier with the metal film, too. 
At this interface, Tg equals T (in the numerics it dif- 
fers marginally since the infinitesimal e is approximated 
by a variable of finite extent). The thickness b enters 
this approximation via a modification of G(r) in Eq. (6), 
which we label G5(r). It still obeys the Laplace equa- 
tion V^G^ = 0, now in the barrier, i.e., for — 6 < z < 0, 
but the Dirichlet condition G(r -^ oo) = is replaced at 
z = —b by Gb{z = —b) = 0. In this way we approximate 
the temperature at the interface between barrier and sub- 
strate by the base temperature of the unbiased sample. 
We refer to this type of approximation as Approx. II in 
the following section. The calculation of G5(r) on a rect- 
angular grid used in the finite- difference approximation 
is explained in the appendix. 



IV. RESULTS 

The model developed in the last section can be used to 
optimize structure layouts intended for the formation of 
narrow break junctions at well defined locations by de- 
liberate electromigation. In this section we present some 
examples. 
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Figure 3. Maps of normalized temperature, T/l'^Ru^ of metallic structure [(a) and (c)], substrate-barrier interface [(b) and 
(d)], and interface between barrier and metallic structure [(e) and (f)]. Colors are according to the scale in units of K/mW. 
Results for three barrier thicknesses are shown: 6=1 nm [(a) and (b)], b = 100 nm [(c), (d), and (e)], and b = 1 ]im (f). (a) to 
(d) are calculated with Approx. I, while for (e) and (f) Approx. II is used. 



First we want to address the question why metallic 
structures on a thicker silicon oxide layer are more likely 
to melt during the electromigration process than sam- 
ples on native Si02. For that purpose we calculate the 
excess-temperature map for the sample geometry shown 
in Fig. 2 for varying barrier height h. For definiteness 
we consider gold structures {km ^ 316W/Km^^) with a 
layer thickness of <i = 30 nm, approximate the thermal 
conductivity of the substrate by the typical value of sil- 
icon at room temperature, kg ~ 150W/Km^^, and vary 
h between 1 nm and 2.5 pm. For the thermal conductiv- 
ity of the barrier we use that of thin Si02 films at room 
temperature, A:^ ^ l.lW/Km^^. 

Figure 3 shows the results for 6 = 1 nm, h = 100 nm, 
and 6=1 pm. Obviously, the temperatures of substrate 
and metallic structure follow each other closely for very 
thin barriers [see Fig. 3 (a) and (b)]. For thicker barriers 
the asymptotic temperature T^ in the leads can easily 
reach values that exceed the maximum temperature for 
thin barriers [compare Fig. 3 (a) and (c), calculated with 
Approx. I, i. e., a 3d heat relaxation in the substrate and 
a vertical laminar flow in the barrier]. The temperature 
is maximal at the center of the wire as expected. A typ- 
ical value for the sheet resistance is R^ = 2 r^, so that 
PRu = 1 mW at around / = 20 mA. For this current 
a moderate temperature rise of T^ax ^ 50 K results at 
the center of the narrow wire, if the barrier thickness 



6=1 nm. For 6 = 100 nm the temperature rises already 
by Tmax ^ 510 K. The substrate temperature below the 
leads, however, does not change significantly for thicker 
barriers due to the fact that the substrate underneath 
long leads is heated by the power density PRu/w'^ in- 
dependent of barrier thickness [compare Fig. 3 (b) and 
(d)]. For i?n = 2 r^ we get a marginal warming by about 
1.6K at / = 20 mA. Furthermore, substrate heating be- 
low the center of the wire is reduced with increasing 6. 
Between 6 = 1 nm and 6 = 100 nm the maximum tem- 
perature rise in the substrate drops from about 30 K to 

10 K. This is a rather low temperature increase as com- 
pared to the rise in the metallic structure, i.e., Approx. 

11 is a better choice for 6 > 100 nm. Figure 3 (e) displays 
the temperature profile in the plain defined by the in- 
terface between the barrier and metal film calculated in 
the framework of Approx. II, i. e. evaluating the 3d heat 
relaxation in the barrier, while the substrate is assumed 
to remain cold. Approximations I and II lead to similar 
results at 6 = 100 nm as can be judged by inspecting Fig. 
3 (c) and (e). Barrier heating is restricted mainly to the 
regions underneath the metallic structure for 6 = 100 nm 
with a small halo only. At larger 6 [see Fig. 3 (f)], how- 
ever, the area of the barrier with a significant tempera- 
ture rise has a considerably larger extent. Nevertheless, a 
pronounced maximum of the temperature in the central 
spot of the structure remains the dominant feature. ^^ 
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Figure 4. Excess temperature T along the central axis of the 
same structure as used in Fig. 3. Note the logarithmic vertical 
scale. A picture of the structure has been underlayed to ease 
matching regions of different current densities with the shown 
curves. L — 1 pm is defined as the length of the narrow wire 
in the center. The thickness h of the Si02 barrier is varied 
between 1 nm and 2.56 pm, the temperature is normalized by 
PRu- Solid lines correspond to Approx. I, while dashed lines 
are calculated using Approx. II. 



Data for further values of h are summarized in Fig. 
4, where the temperature profile along the central axis 
is displayed. The successive rise of the asymptotic lead 
temperature T^ with h for \x/L\ ^ 4 is the most promi- 
nent feature in this figure. Furthermore, one can ob- 
serve that the ratio between the maximum temperature 
in the center and T^ is strongly reduced with increasing h. 
For h = 100 nm results of Approx. I and II are included. 
Both approximations yield similar results at this thick- 
ness, with slightly lower values of T for Approx. II due to 
the more efficient cooling resulting from the temperature 
halo around the metallic structure visible in Fig. 3 (e). 

In a well accepted approximation^ , the amount of gold 
pushed in the anode direction is proportional to both 
the current density and the self-diffusivity of the metallic 
cores. The latter has a strong temperature dependence 
described over a wide range by an Arrhenius law^^ (for 
self-diffusion data on Au see, e. g., Refs. 22-25). No mat- 
ter whether electromgration is dominated by diffusion in 
the bulk, along grain boundaries, or at the surface, a 
sizeable effect sufficient for the deliberate formation of 
voids requires elevated temperatures. As a rule of thumb, 
T should be at least half the melting temperature T^, 
in which case electromigration along grain boundaries is 
the dominant atomic flux mechanism. The effect of elec- 
tromigration is much stronger just below melting and 
is then in general dominated by flux within the bulk. 
However, for a controlled process of gap formation the 
melting point must not be exceeded. In Fig. 5 (a) we 
plot a = \/l/y^^ as a function of 6, where ?/max are the 
maximal values of the curves in Fig. 4 at x = 0, where 
the temperature in the metallic structure is highest. If we 
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Figure 5. (a) Plot of a = lyjRu/To in units of mky/OjK 
as a function of the barrier thickness h. Here / is the current 
needed to heat the structure of Fig. 3 in its center to a tem- 
perature To. (b) Ratio 7 of the first term in Eq. (1) to the neg- 
ative of the term on the right hand side of the same equation 
in the center of the structure, 7 = (dbkm / kb)V'^T / (Ts — T). 
This compares the amount of heat flow in the metal to the 
effect of cooling by the substrate. At large b the structure 
is thermally decoupled in its center, (c) Ratio /3 = Ti/Tmax 
of the temperature Ti at the crossover from the narrow wire 
to the funnel-shaped parts of the leads and the temperature 
Tmax in the center as a function of b. Red dots in (a), (b), 
and (c) are calculated with Approx. I, blue dots with Approx. 
II. (d) Temperature profile along the central axis of two dif- 
ferent structures at 6 = 10, 100 and 1000 nm. The bluish 
and reddish underlay indicate current and voltage probes, re- 
spectively. Red lines correspond to results for the four-point 
configuration shown in Fig. 3 (complete underlay). Blue lines 
are calculated for a two point configuration, i.e., the voltage 
probes are missing (bluish underlay only). For b — 10 nm 
Approx. I and for larger b Approx. II has been used. 



fix the temperature at a value Tq well below the melting 
point (let's say at Tq ~ 1000 K for gold) /q = ol\/Tq/ R^ 
directly gives a safe current limit for the specific struc- 
ture. This safe current limit is considerably reduced when 
b is increased as shown in 5 (a). For a native oxide barrier 
of 6 ~ 1 nm and Ru ~ 2 ^2 a current of almost 100 mA 
can be tolerated. For the thermally grown oxide layer of 
our experiments [b = 400 nm) the safe current limit is by 
a factor of four smaller. This makes controlled electromi- 
gration considerably more difficult. The feedback loop as 
implemented by us and others^^ uses a resistance change 
to judge the appropriate current density. This procedure 
however might lead to current densities that exceed the 
limit set by the melting point of gold. 

It is instructive to compare the size of the first term in 
Eq. (1), i. e., the heat transport in the metallic structure, 
to the size of the right-hand side of the same equation, 
i. e., the heat flow to the substrate. This is done in Fig. 5 
(b) for the central point x = of the structure where the 



temperature for a given current / is highest. The ratio 

7=(d6fc„/fcb)(V2T/(T,-r)) 

of these two terms is considerably smaher than 1 for 
b < 5nm, i.e., the center of the structure is efficiently 
cooled via the substrate. For larger 6, however, 7 rises 
quickly indicating that the coupling to the substrate be- 
comes increasingly inefficient and cooling has to rely on 
the heat flow along the metallic leads. The barrier thick- 
ness at which this decoupling sets in is expected to be 
reduced with decreasing width of the notch in the center. 
Thus, if void formation by electromigration occurs and 
the current path through a narrowing constriction, tem- 
perature is dominated more and more by the heat flow 
along the structure. This leads to an even stronger ten- 
dency towards melting which has to be counteracted by 
increasingly careful settings in the feedback-control loop 
and necessitates good thermal coupling from the begin- 
ning. In summary, our analysis gives clear evidence why 
controlled electromigation on thick oxide barriers is more 
likely to fail. 

Figure 5 (c) displays the ratio (3 of the temperatures 
at the wire end (i.e., at the crossover from the narrow 
region to the funnel-shaped parts of the leads) and at 
the center, as a function of b. Ideally, the temperature 
peaks so strongly at the center that electromigration acts 
only there with reasonable strength. This is difficult to 
achieve, though. Figure 5 (b) indicates that an optimal 
barrier thickness 6opt exists at which f3{bopt) ~ 0.5. This 
factor-of-two reduction of the temperature is not suffi- 
cient as demonstrated by the experimental result pre- 
sented in Fig. 1 (a). As indicated in the inset of Fig. 2 
the current density peaks not only at the notch in the 
center but at concave corners as well. Electromigration 
leads to material flux in regions of high current density. 
The formation of voids or hillocks, in turn, requires some 
kind of inhomogeneity that causes a divergence of ma- 
terial flux. Such an inhomogeneity might be simply due 
to the statistical distribution of grains comprising the 
metallic structure, but regions of high temperature gra- 
dient lead for sure to a divergence in atomic flux as the 
delivery of material from cold regions is less efficient than 
delivery of material to the hotter parts. In the structure 
studied here the temperature gradient is large at the wire 
end and vanishes in the center where the notch is located. 
In the latter position material transport tends to be bal- 
anced and thus the wire cross-section does not change to 
first order. Only at the corner at the cathode end void 
formation is favored. 

Figure 5 (d) compares calculations for two different 
structures, namely the four-probe configuration with 
source and drain contacts and two voltage probes stud- 
ied so far, and a two-probe configuration with no voltage 
probes. Clearly the temperature is reduced in the vicin- 
ity of the current-free voltage probes. However, their 
cooling effect is confined to a small region whose extent 
is set by the relaxation parameter A = yWkjrJkiy in Eq. 
4. The design of cooling fins would require to place them 
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Figure 6. (a) Equipotential lines for a ring structure. The 
voltage drop between consecutive lines is /SV = 0.1 i?n • I 
and the total voltage drop between the outermost lines is 
/SV — 13.8i?n •^- The enlarged detail shown in the ring center 
displays the power density in units of I^Ru/w^ (see color 
scale on the right bottom), (b) Temperature profile along 
the blue line in (a) for various barrier thicknesses. Solid lines 
are calculated with Approx. I and dashed lines with Approx. 
II. (c) Temperature map of the ring structure with a barrier 
thickness oib — 1 nm according to Approx. I. 



within a distance set by A. For the model parameter used 
in our calculations A r^ 290 nm, 930 nm, and 2.9 pm for 
6 = 10 nm, 100 nm, and 1 pm, respectively. The voltage 
probes in our test structure are within the range of the 
central notch for 6 = 1 pm, only. Consequently, the tem- 
perature is slightly reduced in the center in this case as 
can be seen in Fig. 5 (d). Note, however, that the tem- 
perature in the center of the structure is peaking on top 
of the asymptotic temperature T^ and this temperature 
scales quadratically with 1/w^ the width of the source 
and drain leads. Thus, a stronger reduction can be ob- 
tained by making the leads as wide as possible. 

To demonstrate the application of our numerical sim- 
ulations for the optimization of the sample layout with 
respect to control over the location where gap formation 
sets in, we design a structure in which the current path 
branches into a loop (see Fig. 6). The final goal with 
this specific structure is to produce two gaps, one in each 
loop arm. If bridged by molecules, such a sample could 
be used to test for the quantum coherence of electron 
transport through molecules by measuring the magneto- 
conductance. Coherence would lead to the appearance 
of periodic oscillation of the conductance in response to 
the magnetic flux penetrating the area between the loop 
arms^^. In the design of Fig. 6 (a) we try to avoid unin- 
tentional current density enhancements by rounding all 
concave corners. A step in the width of the current path 
is introduced only at two positions. The inset in the ring 



center of Fig. 6 (a) demonstrates that this leads to a con- 
siderable current-density enhancement at singular points. 
At the same time we would like to have a large temper- 
ature gradient at this location. To this end we make the 
current path twice as narrow on the anode side of the 
ring structure and thus the power density by a factor of 
four larger. The investigation of coherent transport re- 
quires a high aspect ratio of the open area between the 
loop arms and the part of the ring area covered by metal. 
For this reason, the step in the ring width is placed close 
to the cathode side of the structure. Figure 6 (c) shows 
the result of Approx. I for a thin barrier (6 = 1 nm). The 
temperature changes abruptly at the transition from the 
wide to the narrow part of the ring thus making the con- 
cave corner at this position a preferred starting point for 
void formation. In Fig. 6 (b) we display the temperature 
profile along the path drawn as a blue line in Fig. 6 (a) for 
various barrier thicknesses. The temperature gradient is 
large at point B for a thin barrier as intended. The tem- 
perature reaches to first approximation almost constant 
values to the left of point A and between point C and D 
for 6 = 1 nm on account of a short relaxation parameter 
A ~ 290 nm. With increasing thickness 6, however, Joule 
heating leads to a more evenly distributed temperature 
due to the rise in A ex \/h. As a result, the temperature 
gradient at point B gets more and more shallow mak- 
ing the tendency to void formation at the corner weaker. 
Still, electromigration is boosted by the current-density 
rise at this location which is thus distinguished from the 
rest of the structure. Nevertheless, the statistical nature 
of grain boundaries easily marks different locations along 
the line from B to D where void formation might start. 
Hence, it is not clear from the outset that the corner 
close to point B will win for larger h. This again stresses 
that reliable void formation is favored by good thermal 
coupling via a thin barrier to a substrate of high thermal 
conductance. 



perature is governed only by heat flow along the leads 
connecting the microstructures to regions of lower cur- 
rent densities. Thus, only thin barriers allow the large 
current densities required for void formation on a reason- 
able time scale without melting large parts of the struc- 
ture by Joule heating. Therefore, good thermal coupling 
is an important prerequisite for well controlled electromi- 
gration. On the other hand, thin barriers are not always 
an option in practice. The native oxides on silicon wafers 
have been reported to be as thin as 1 nm. While such a 
thin barrier yields a good thermal contact it bears the risk 
of undesirable electrical leakage currents through the sub- 
strate. For a reliable electrical insulation of the structure 
a barrier of several 100 nm thickness is mandatory. In our 
model calculation we considered cooling the structure by 
coupling it to a cold substrate. If electromigration has 
to be performed in vacuum for some reason, this might 
be the only option. In this case a thermally well con- 
ducting insulator is required. Another option would be 
to operate electromigration in a controlled atmosphere. 
We have, e.g., performed electromigration experiments 
in the gas phase above the liquid level of a Helium de- 
war. The structure heated by a high current density 
leads to convection of the cold gas which in turn gives 
rise to an efficient heat release. The yield of sufficiently 
small nanogaps considerably increased in this case even 
on thick Si02 barriers. Whether the resulting gaps are 
clean and not contaminated with metallic nano-clusters 
is an open question at this moment^^. 

Besides efficient heat release, a good layout design is 
needed to predefine the position of void formation. The 
latter starts preferably at positions where high current 
densities meet large temperature gradients. It is thus 
advisable not to place the intended gap position in high 
symmetry points of the layout where temperature gra- 
dients are expected to vanish but to mark it with sharp 
corners which can be utilized for a considerable current 
density enhancement. 



CONCLUSION 
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Appendix: Calculation of G'b(r) 

We calculate ^^(r) by discrete cosine transform (DOT) 
on a cube of A/" x TV x TV^ cells (N=1024 turns out to be 
sufficient) with the basic approach 



-1 ^z 



Gb{x,y,z)= Y, Yl Gb{X,Y,Z)cos{7Tx{X^l/2)/n)cos{7Ty{Y^l/2)/n)cos{7T{z^l/2){Z^l/2)/N,), (A.l) 



x,y=0 z=0 



where n = N/2 and the arguments of the co- 
sine functions are chosen to assure Gb{^n^y^z) = 
Gb{x,±n,z) = Gb{x,y,N, - 1/2) = 0, Gb{l,y,z) = 
Gb{-1, y, z), Gb{x, 1, z) = Gb{x, -1, z), and Gb{x, y, 0) = 
Gi,{x^y^—1). Together with the discrete version of the 
Laplace equation, 



D'^G = 5r 



(A.2) 



this describes a single source of unit strength at {x = 
0,^ = 0,z = —1/2). D'^ is the discrete version of the 



Laplace operator V , 6i 



000 



1, and J^ 



for all other 



combinations of the indices. From Eq. (A.2) one gets 



Gb{X,Y,Z) 



q{X,Y,Z) 



4(cos(7r(X + l/2)/n + cos(7r(Y' + l/2)/n) + cos{tt{2Z + 1)/N,)) 



(A.3) 



where g(X, F, Z) is the DCT of 6xyz defined with the 
same arguments in the cosine functions as in Eq. (A.l). 
Thus the calculation of G5(r) involves two DCT opera- 
tions and normalization of the intermediate result accord- 
ing to Eq. (A.3). The index b enters into the calculation 



I 

via the definition of N^ = b/h, where h is the discretiza- 
tion size in the finite- difference method. For Approx. I we 
need to know G(r) = lim5^oo ^^(r). We approximate it 
by setting G(r) ^ G512, which is sufficiently accurate due 
to the asymptotic 1/r nature of the analytical solution. 
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